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Abstract. We compare the hot-star wind models calculated by assuming older solar-abundance determination with models 
calculated using the recently published values derived from 3D hydrodynamical model atmospheres. We show that the use of 
new abundances with lower metallicity improves the agreement between wind observation and theory in several aspects. (1) 
The predicted wind mass-loss rates are lower by a factor of 0.76. This leads to better agreement with mass-loss rates derived 
from observational analysis that takes the clumping into account. (2) As a result of the lowering of mass-loss rates, there is 
better agreement between the predicted modified wind momentum-luminosity relationship and that derived from observational 
analysis that takes the clumping into account. (3) Both the lower mass fraction of heavier elements and lower mass-loss rates 
lead to a decrease in opacity in the X-ray region. This influences the prediction of the X-ray line profile shapes. (4) There is 
better agreement between predicted P v ionization fractions and those derived from observations. 
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1. Introduction 

Massive hot stars lose a substantial part of their mass via their 
stellar winds. Wind parameters, especially the wind mass-loss 
rate (the amount of mass lost by the star per time unit), are the 
functions of basic stellar parameters, for example the stellar lu- 
minosity (see, e.g., Kudritzki & Puis 2000; Krticka & Kubat 
2006, for a review). The wind mass-loss rate is the most impor- 
tant wind parameter, because the presence of strong hot-star 
winds influences not only the emergent spectra, but also, for 
example, stellar evolution (e.g., Chiosi & Maeder 1986). 

Since hot-star winds are mainly accelerated by the absorp- 
tion (or scattering) of radiation in lines of heavier elements 
(like carbon, nitrogen, or iron), the radiative force and wind 
mass-loss rate also depend significantly on metallicity. With 
increasing metallicity, the mass-loss rate and the wind ter- 
minal velocity increase (Kudritzki & Puis 2000; Vink et al. 
2001; Krticka 2006). Frequently, the solar abundance values 
of Anders & Grevesse (1989) are used for the wind-model pre- 
dictions. However, the solar metallicity was recently reduced 
following the calculations by Asplund et al. (2005) based on 
improved (self-consistent) treatment of the outer solar con- 
vective region (using 3D hydrodynamical model atmospheres, 
Nordlund & Dravins 1990; Steffen et al. 1995; Asplund 2004). 
The difference between older abundance determinations and 
recent, more realistic ones is very striking, as the estimate of 
the mass fraction of heavier elements decreased from 0.0194 
to 0.0122. Consequently, it would be interesting to test the in- 
fluence of the new solar abundance determination on the wind 
models. To do so, we compare wind models calculated with 
both old and new abundance determinations. 



2. Basic assumptions 

2.1. Wind models 

We used our own NLTE wind models to predict the hot- 
star wind structure (Krticka & Kubat 2004, recent improve- 
ments in our calculations, especially the inclusion of the Auger 
ionization, will be described elsewhere). Our models enabled 
us to self-consistently solve the hydrodynamical equations in 
the spherically-symmetric, stationary, radiatively-driven stel- 
lar wind. The radiative force (in the Sobolev approximation, 
Castor 1974) and the radiative heating/cooling term were cal- 
culated using the solution of NLTE equations, for which we 
mainly used the atomic data available from the Opacity and 
Iron Projects (Seaton et al. 1992; Hummer et al. 1993). This 
enabled us to change the abundance of each element separately 
according to its new value, and not only via the common pa- 
rameter called the metallicity. Using our method we were able 
to predict wind density and velocity structure, namely the wind 
mass-loss rate M and the wind terminal velocity iw 

2.2. Model stars 

For our analysis we selected O stars with T e ff < 40 000 K (see 
Table 1), for which the stellar parameters were derived using 
the models with line blanketing. Effective temperatures and 
radii were taken from Repolust, Puis, & Herrero (2004, here- 
after R04), Markova et al. (2004, hereafter M04), and Martins 
et al. (2005, hereafter M05). Stellar masses were obtained us- 
ing evolutionary tracks either by us (using Schaller et al. 1992 
tracks) or by M05. 
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Table 1. Stellar parameters of the O stars. 
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Data sources: M04 - Markova et al. (2004), R04 - Repolust et al. 
(2004), M05 - Martins et al. (2005) 



Table 2. Comparison of wind parameters calculated assuming Anders 
& Grevesse (1989) and Asplund et al. (2005) solar abundances. 
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3. Wind models with Asplund, Grevesse, & Sauval 
solar abundance determination 

With decreasing metallicity the radiative force decreases 
and, consequently, the wind mass-loss rate also decreases. 
Moreover, the decrease in metallicity leads to a slight lower- 
ing of the wind terminal velocity (Kudritzki & Puis 2000; Vink 
et al. 2001; Krticka 2006). Both these effects are apparent in 
Table 2, where we compare wind parameters calculated using 
different determinations of solar abundances (see also Fig. 1). 
On average, the mass-loss rate changes by a factor 0.76. This 
roughly corresponds (the abundance changes were not uniform 
for all elements) to the dependence M ~ Z° 60 , where Z is the 
mass fraction of heavier elements. This dependence is slightly 
lower than the one reported by Vink et al. (2001) and Krticka 
(2006). Note, however, that these results are not directly com- 
parable, since the relative change in the abundance of differ- 
ent elements is not uniform in this study, whereas results by 
both Vink et al. (2001) and Krticka (2006) were obtained for 
a uniform relative change in the abundance. On average, the 
wind terminal velocity also decreases with decreasing metallic- 
ity (see Fig. 1). However, note that the change in the terminal 
velocity due to the decrease in the metallicity is mostly lower 
than the error in the terminal velocity determination. 

4. Implications of new abundance determinations 

The next logical step would be to compare our predicted wind 
parameters (especially the mass-loss rates) with the quantities 



models with new solar abundances 
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Fig. 1. Comparison of the ratio of the terminal velocity to the escape 
velocity calculated using terminal velocities either from theoretical 
calculations (with different sets of abundance determinations) or de- 
rived from observations (values adopted from Puis et al. 1996, M04, 
M05). LSL denotes the average ratio derived by Lamers et al. (1995). 

derived from observation. However, such a comparison is com- 
plicated by the effect of clumping. 

4.1. Wind clumping 

There are two ways to derive wind mass-loss rates. Since many 
spectral features are sensitive to wind density (or mass-loss 
rate), it is possible to derive their values from observations us- 
ing wind models where the mass-loss rate is a free parameter. 
These values are sometimes called 'observed', but we prefer 
to call them 'derived from observations'. On the other hand, 
the mass-loss rate can be calculated ab initio from dedicated 
wind models where stellar parameters (the effective tempera- 
ture, mass, radius, and metallicity) serve as input parameters. 
This value of the mass-loss rate is usually called 'predicted' . It 
is clear that both ways of determining the wind mass-loss rate 
should be in close agreement. 

However, recent spectroscopic studies (Bouret et al. 2003; 
Martins et al. 2005; Puis et al. 2006) indicate that there may 
be significant disagreement between the wind mass-loss rates 
derived from observations and predicted ones. The problem is 
that the emergent spectrum is, in many cases, not directly sen- 
sitive to the wind mass-loss rate M, but to the wind density 
squared. This means that, if it happens that the wind is struc- 
tured on small scales ("clumped"), then the lines originating in 
the clumped wind with a relatively low mean wind density may 
mimic those with a much higher mass-loss rate. If the winds 
are really significantly clumped in all regions from which we 
observe the emergent radiation, as seems to follow from the 
comparison of synthetic wind spectra with the observed ones 
(e.g., Bouret et al. 2003; Martins et al. 2005), then the real wind 
mass-loss rates are significantly lower than the predicted ones. 

From the studies of wind spectra that include clumping 
(M05), it follows that wind mass-loss rates are roughly lower 
by a factor ~ 0.3 than those derived from the standard wind 
models of Vink et al. (2001). Consequently, the reduction of 
wind mass-loss rate by a factor 0.76 due to the reduced metal- 
licity represents an improvement in the agreement between 
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Fig. 2. Comparison of modified wind-momentum calculated for in- 
dividual stars in this paper assuming solar abundances by Anders & 
Grevesse (1989) with more recent lower ones by Asplund et al. (2005). 
We also plot the linear fits of the theoretical predictions, derived with 
older abundances (Krticka & Kubat 2004; Vink et al. 2001), and values 
derived from observations both by M05 (individual stars for which the 
clumping was taken into an account) and by R04 (linear fit for giant 
and dwarfs, for which the clumping was not taken into an account). 

wind mass-loss rates derived from observation (that take the 
clumping into account) and predicted ones. 

4.2. Modified wind momentum-luminosity relationship 

It follows from both wind theory and observations that the 
value of the modified wind momentum Mv m (R/Rq) 1/2 de- 
pends mainly on the stellar luminosity (Kudritzki & Puis 2000, 
and references therein). Especially the mass-dependence nearly 
cancels out. Consequently, the wind momentum-luminosity re- 
lationship is suitable for comparing observations and theory. 
Here we compare this relationship derived using theoretical 
models with both older and new abundance determinations 
with values from M05, who derived wind mass-loss rates from 
the observed data and who account for clumping, and R04, who 
neglected the possibility of wind clumping (see Fig. 2). 

The predicted wind momentum is lower for models calcu- 
lated using improved solar metallicity determination, mainly 
due to the decrease in wind mass-loss rate with metallicity 
and, to a lesser extent, due to the lower terminal velocities 
that are derived. The theoretical modified wind momentum- 
luminosity relationship calculated using improved abundances 
has the form 

log[M Vca (R/RQ) l/2 ] = 1.91 log(L/L Q ) + 18.1 (CGS). (1) 

Our results to some extent improve the agreement between the 
theory and observational analysis that takes the clumping into 
account (M05 and also Fig. 2). However, they are not able to re- 
move the difference between the theory and observations com- 
pletely. Note that there is good agreement between theoreti- 
cal predictions derived using Anders & Grevesse (1989) abun- 
dances (Vink et al. 2001; Krticka & Kubat 2004, and the empty 
circles in the Fig. 2) and observations by R04. Consequently, 
if the influence of clumping on the observed wind spectra is 
negligible, then the wind models with improved abundance de- 



Table 3. Comparison of continuum opacity (opacity per unit of mass 
in units of cm 2 g _1 averaged for the radius 1.5 Rq < r < 5R Q ) in the 
X-ray region for several different wavelengths calculated for the two 
abundance determinations (for HD 24912). 



wavelength [A] 8.1 12.0 15.8 20.0 24.0 
Anders & Grevesse (1989) 31 83 141 228 154 
Asplund et al. (2005) 20 52 91 150 125 



termination predict lower wind momentum than those derived 
from observations. 

4.3. X-ray continuum optical depth 

Another painful problem of the hot-star wind theory is con- 
nected with the shape of X-ray line profiles. While the wind 
theory based on currently available theoretical mass-loss rates 
predicts asymmetric X-ray line profiles due to the continuum 
absorption within the wind (see, e.g., Owocki & Cohen 2006; 
Oskinova et al. 2006), the observations mostly show symmetric 
X-ray line profiles. This indicates that the wind's optical depth 
in continuum for X-ray frequencies is small in the region where 
X-rays form. 

To infer the impact of new solar abundances on the X-ray 
line profiles, we compared the mass-absorption coefficients cal- 
culated with solar abundances from both Anders & Grevesse 
(1989) and Asplund et al. (2005). Apparently, since heavier el- 
ements dominate the continuum opacity in the X-ray region 
(especially due to Auger transitions), the mass absorption coef- 
ficient decreases with decreasing metallicity (see Table 3). The 
decrease in the optical depth is not equal for all wavelengths 
because different elements dominate the opacity for different 
wavelengths and the relative change of the metallicity is not 
the same for all elements. The influence of new solar abun- 
dance determination on the total continuum opacity in the X- 
ray region is amplified by the fact that the wind density also 
decreases with decreasing metallicity. 

As a result of these effects, the X-ray optical depth is lower 
roughly by a factor of 2. The use of lower X-ray optical depth 
in the calculation of Owocki & Cohen (2006) leads to better 
agreement between predicted and observed X-ray line shapes. 
Conversely, the effect of scattering in these lines (Ignace & 
Gayley 2002) might decrease with metallicity. 

4.4. Phosphorus line profiles 

The ion P v is assumed to be a dominant phosphorus ion in 
the stellar wind. Since its abundance is very low, its resonance 
lines as observed in the hot star wind spectra are generally un- 
saturated. Consequently, these lines are assumed not to be in- 
fluenced by clumping and are ideal for determining the wind 
mass-loss rate. This was done by Fullerton et al. (2006), who 
showed that the mass-loss rate determined from the P v lines 
is much lower than when derived from the Ha line or radio 
emission. 

This discrepancy can have two explanations: either Pv is 
not a dominant phosphorus ion or the winds are clumped and 
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their mass-loss rates are much lower. However, the phospho- 
rus line profiles are sensitive to phosphorus abundance and to 
the overall stellar metallicity (due to the dependence of wind 
mass-loss rates on the metallicity). Both these effects lead to 
the weakening of P v line profiles. The phosphorus abundance 
changes by a factor of 0.8 1 when using a new abundance deter- 
mination and, together with the decrease in wind mass-loss rate 
caused by the decrease in the metallicity, leads to the change in 
Pv line optical depth by a factor of about 0.62. Although this 
shift does not explain the main difference between theory and 
observations (possibly also the influence of clumping on the 
ionization balance has to be accounted for, Puis et al. 2006a), 
the use of new abundance determinations shifts the theoretical 
predictions closer to the observational results. 

5. Discussion and conclusions 

We have compared the wind models calculated using the older 
solar abundance determination (Anders & Grevesse 1989) and 
the solar abundances derived using 3D hydrodynamical model 
atmospheres of Asplund et al. (2005). We have shown that the 
use of new abundances (which means lowering in the metallic- 
ity, but not uniformly for all elements) improves the agreement 
between observation and theory in several different aspects: 

- Due to lower metallicity the predicted wind mass-loss rates 
are lower, consequently new models agree better with the 
mass-loss rates derived from observational analysis with 
clumping taken into an account. 

- As a result of the lowering of mass-loss rates and slight low- 
ering of terminal velocities, there is better agreement be- 
tween the predicted modified wind momentum-luminosity 
relationship and that derived from observational analysis 
with clumping taken into account. On the other hand, if 
the clumping does not significantly influence the observed 
spectra, then the original good agreement between pre- 
dicted wind momentum and the one derived from observa- 
tion is worsened with models using improved solar abun- 
dance determination. 

- Both the lower mass fraction of heavier elements and lower 
mass-loss rates lead to a decrease in the opacity in the X- 
ray region. This change is the most significant one because 
it basically combines the influence of two different effects 
(when using predicted mass-loss rates). 

- Since both phosphorus abundance and the mass-loss rate 
are lower, there is better agreement between the predicted 
Pv ionization fraction and the one derived from observa- 
tions. However, lowering the wind density may partly com- 
pensate for this effect due to modification of the ionization 
structure. 

We have to keep in mind that the actual metallicity of Galactic 
O stars may be different from the solar one. Moreover, it may 
be modified during the stellar evolution due to the rotational 
mixing (e.g., Maeder & Meynet 2001), so our results have to 
be understood as a manifestation of the consequences of abun- 
dance changes in the theoretical wind models. 

We can conclude that, although there is still not complete 
agreement between theory and observations, the use of more 



realistic solar abundance determinations may be a partial way 
to solve the problem because it improves the agreement in all 
aspects studied here. 
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